Modified f(R) gravity from scalar-tensor theory and inhomogeneous EoS dark energy 
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The reconstruction of f(R)-gravity is showed by using an auxiliary scalar field in the context of 
cosmological evolution, this development provide a way of reconstruct the form of the function f(R) 
for a given evolution of the Hubble parameter. In analogy, f(R)-gravity may be expressed by a perfect 
fluid with an inhomogeneous equation of state (EoS) that depends on the Hubble parameter and 
its derivatives. This mathematical equivalence that may confuse about the origin of the mechanism 
that produces the current acceleration, and possibly the whole evolution of the Hubble parameter, 
is shown here. 
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I. INTRODUCTION 



Ever since the accelerated nature of the dynamics of the Universe was discovered in 1998 by the observations of 
SN la, a lot of theoretical descriptions have been proposed such that the observational data is satisfied. They have 
established that the EoS parameter w for the fluid that governs the Universe is close to -1. The majority of these 
models are included in the so-called dark energy models, where the origin of the dark energy is a scalar field or a perfect 
__■ fluid with an inhomogeneous equation of state(EoS), which should be the responsable for the current acceleration, and 
\ may be even on the whole expansion history of the Universe. On the other hand, the modified f(R) theories of gravity 
"^T* ' avoid the need to introduce dark energy, and may give an explanation about the origin of the current accelerated 
Oh' expansion and even on the expansion history of the Universe(for recent reviews, see [l| and [2(]). 

In this sense, the cosmic acceleration and the cosmological properties of metric formulation f(R) theories have been 
studied in Refs.[^|-(33|. Recently the main focus has been improviding a f(R)-theory that reproduces the whole history 
of the Universe, including the early accelerated epoch (inflation) and the late-time acceleration at the current epoch 
(see [H]|-|38j|), where the possible future singularities have been studied in the context of f(R)-gravity(see [38j|). It is 
important to remark that the main problem that this kind of theories found at the begining of its development was 
the local gravitational test; nowadays several viable models have been proposed, which pass the solar system tests 
\ and reproduce the cosmological history (see [35l|-(47j). 
I . In the present paper, the reconstruction of f(R)-gravity is shown, to be possible in the cosmological context by using 
an auxiliary scalar field and then various examples are given where the current accelerated expansion is reproduced 
and also the whole history of the Universe. This kind of reconstruction is well done by using an scalar field without 
kinetic term, which differs of the kind of quintessence reconstruction models (see [48| and [49|]) where the scalar field 
presented has a non-zero kinetic term. Also it is investigated the analogy between the so-called dark fluids, whose EoS 
is inhomogeneous and which have been proposed as candidate of dark energy (see [I(|-[55|]), and the f(R)-theories, a 
reconstruction of such kind of theories is shown by using such type of perfect fluids. 

The paper is organized in two sections, in the first one the reconstruction of f(R)-gravity is obtained by using scalar- 
tensor theory. In the second section, the analogy between the considerations of dark fluids with an EoS depending on 
the Hubble parameter and its derivatives, and f(R)-gravity is considered. 



II. RECONSTRUCTION OF F(R)-GRAVITY 

In this section, it will be shown how f(R) theory may be reconstructed in such a way that cosmological solutions 
can be obtained. Let us start with the action for f(R)-gravity: 
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S= / d*xV=g(f(R)+L m ) . (1) 
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Here L m denotes the lagrangian of some kind of matter. The field equations are obtained by varying the action on 
<7 M „, then they are given by: 

R^f'(R) - \g^f{R) + g^of(R) - v^ v f'{R) = yT<?> , (2) 

where is the energy-momentum tensor for the matter that filled the Universe. We assume a flat FRW metric: 

3 

ds 2 = -dt 2 + a 2 (t)^dx 2 . (3) 



i=l 



Then, if is a perfect fluid, the modified Friedmann equations for the Hubble parameter H(t) — -, take the form: 

l -f{R) - 3(H 2 + H)f(R) + mf"(R)(H 2 H + HH) = ^p m , 

\f{R) (3H 2 + H)f'(R) - Df(R) = ~^ Pm , (4) 

where the Ricci scalar is given by R = 6(2H 2 + H). Hence, by the equations (H|), any cosmology may be reproduced 
for a given function f(R). Nevertheless, in general it is very difficult to get an exact cosmological solution directly 
from (j4}. It is a very useful technique developed in [B| and [34j |. where an auxiliary scalar field without kinetic term 
is introduced, then the action Jl} is rewritten as follows: 

s = J dW=sTO)#+QW>)) , (5) 

where the scalar field (f> has no kinetic term. By variation on the metric tensor g^ v the field equation is obtained: 

\g^ (P(ct>)R + Q{4>)) + P{ct>)R»v + 9^UP{4>) - V M V„P(0) = yT^> . (6) 

The action ([5]) gives an additional equation for the scalar field <fi, obtained directly from the action by varying it with 
respect to 4>: 

p'{4>)R + Q'{4>) = , (7) 

here the primes denote derivatives respect <p. This equation may be resolved with the scalar field as a function of R, 
cj) — <f)(R), and then, replacing this result in the action , the action ((T|) is recovered, 

f(R)=P(cf > (R))R + Q(4,(R)) . (8) 

Hence, any cosmological model could be solved by the field equation ([6|), and then by ([7|) and ([8]) the function f(R) 
is obtained. For the metric ([3]), the Friedmann equations read: 

SH ^dT + m2m + l m ~ P ~Y = ' 

2*™ + AH^l + (AH + W 2 )P{c») +P = 0. (9) 
dt z dt 

We redefine the scalar field such that it is chosen to be the time coordinate <f> = t. The perfect fluid define by the 
energy-momentum tensor Tjllf 1 may be seen as a sum of the different components (radiation, cold dark matter,..) 
which filled our Universe and whose equation of state (EoS) is given by p m — w m p m , then by the energy momentum 
conservation p m + 3H(l + w m )p m — 0, it gives: 

p m = Pmo exp f-3(l + w m ) J dtH(t)j . (10) 



Hence, taking into account the equations ([9]) and (|10|) the Hubble parameter may be calculated as a function of the 
scalar field (j>, H = g((f). By combining the equations ©, the function Q(4>) is deleted, and it yields: 



- 2 5 (0)^ + + (1 + w m ) exp 



-2,(1 + w m ) / dcf>g(cf>) 



= 0. (11) 
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By resolving this equation for a given function P(4>), a cosmological solution H(t) is found, and the function Q(tf>) is 
obtained by means of the equation given by ([5]): 

Q{4>) = -6( ff (0) f ptf) - 6g(^)^- ■ (12) 

If we neglect the contribution of matter, then the equation (jTTJ) is a first order differential equation on g(4>), and it 
can be easily resolved. The solution found is the following: 

9 {<t>) = -y/Ptf) J ^ 2 p2/ ( 3^ } + kP ^ > ( 13 ) 

where fc is an integration constant. As an example to show this construction, let us choose the following function 
that, as it is showed below, reproduce late-time acceleration: 

P(0) = (f) a , where a > 1 . (14) 

Then, by the result l|13p . the following solution is found: 

a + 2 

where fc is an integration constant. By the expression (TH?)) , the function Q(</>) is given by: 



Q(4>) = -6 



a(2a + l)\^3^ a 2 (a-l)(2a+l) a 



(fc0) 2 + Ik + v — y ^ + 



a + 2 / r (a + 2) 5 

The function (|15[) gives the following expression for the Hubble parameter: 



(16) 



^) = H^ + ^ili. (17) 
a + 2 t 

This solution may reproduce a Universe that passes through two phases for a conveniently choice of a. For small 
times the Hubble and the acceleration parameter take the expressions: 

Hit) ~ , 
W a + 2 t ' 

a a(a— 1) / a(a — 1)\ 1 

a 



a + 2 V a + 2 J t 2 



(18) 



where if 1 + \/3 < a < 2, the Universe is in a decelerated epoch for small times, which may be interpreted as the 
radiation/matter dominated epochs. When t is large, the Hubble parameter takes the form: 

H(t) = kt a ' 2 . (19) 

This clearly gives an accelerated expansion that coincides with the current expansion that Universe experiences 
nowadays. Finnally, the expression for f(R) © is calculated by means of (fT4")l and (|16p . and by the expression of the 
Ricci scalar R = 6(2g 2 (0) + g((f>)), which is used to get <j>(R). For simplicity, we study the case where a = 2, which 
gives: 



V 2 

By inserting this expression into (jT4j) and (| 16[) . the function f(R) is obtained 



R-2k±^R(T^2k) ^ 



/ s r , , x , N1 R-2k± JR{1 - 2k) , , 

/(i?) = [J2 - 6(*(fc + 1) + 5/2)] ^— i i + const . (21) 

Thus, with this expression for the function f(R), the current cosmic acceleration is reproduced with the solution 
dTTjl . In general, as it is seen in the following example, it is very difficult to reconstruct the function f(R) for the 
whole expansion history, and even more difficult for the kind of models that unify inflation and cosmic acceleration, 
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in this cases it is convenient to study the asymptotic behaviour of the model, and then by resolving the equations, 
the expression for f(R) is obtained . 



As a second example, one could try to reconstruct the whole Universe history, from inflation to cosmic acceleration 
by the f(R)-gravity as is made in (34 ] — [37j| - In this case, we proceed in the inverse way than above, by suggesting a 
function g((f>), and trying to reconstruct the expression f(R) by calculating P(4>) and Q{4>) by means of equations 
©■ We study an example suggested in [56j, where: 



For this function the Hubble parameter takes the form H(t) = + rrr?- To reconstruct the form of f(R) we have 
to resolve the equation pip . For simplicity we study the assimptotically behaviour of H(t) in such a way that allow 
us to resolve easily the equation (jTTJ) for P{4>). Then, for small t (t << t s ), the Hubble and acceleration parameters 
read: 

, s Hi a H x ( H x 2H X 2H Q \ 

H ^~1F a--w[-W--T + t^i) ■ (23) 

As it is observed, for t close to zero, - > 0, so the Universe is in an accelerated epoch during some time, which may be 
interpreted as the inflation epoch, and for t > 1/2 (t << t s ), the Universe enters in a decelerated phase, interpreted 
as the radiation/matter dominated epoch. The equation pip for P{4>), neglecting the matter component, is given by: 

d 2 P(4>) H x dP{4>) 4ffi 

^^~W^'^ m ^°- (24) 

The solution of is: 



_, . 20k 120fe ,_, . 240fc . 120fc 



p(4>) = fcr 4 + -^r 3 + -j^r 2 + -j^r 1 + - (25) 



where k is an integration constant. The function Q(4>) given by HI]), takes the form: 

120fe ,_ 



Q(4>) = -6§ 



kH^-" + 6kcj)- 



(26) 



By using the expression for the Ricci scalar, the relation (f>{R) ~ (12_ff^/i?) 1 / 4 is found, then, the function f(R) for 
small values of <f> is approximately: 

f(R) - ^ ■ (27) 

Hence, by the expression (|27p the early cosmological behaviour of the Universe (|23"]) . where a first accelerated epoch 
(inflation) occurs and after it, a decelerated phase comes (radiation/matter dominated epochs), is reproduced. Let 
us now investigate the large values for t (t close to the Rip time t s ). In this case the Hubble and the acceleration 
parameters for the solution (f2"2"|) take the form: 

9{4>) ~ T r -> ~ 7 7 and . 28 

t s - t s - 1 a [t s - ty 

As it is observed, for large t the solution (|2"2"|) gives an accelerated epoch which enters in a phantom phase (H > 0) 
and ends in a Big Rip singularity at t — t s . In this case the equation for P((j)) reads: 

d 2 P{^) H Q dP{^) 2H 

- —4, — + Mf PW = ° ' (29) 



The possible solutions of the equation (|29|) depend on the value of the constant Ho, as it follows: 
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1. If Hq > 5 + 2\/6 or Ho < 5 — 2^/6, then the following solution for P((j>) is found: 

P(<f>) = A(t a - 4*)°+ + B(t. - <f>) a - , 



where a± 



H Q + l± ^0(^0 + 10) + ! 



(30) 



Then, through the expression (fl"2"|) , the function Q(0) is calculated. In this case, for t close to t s , the Ricci scalar 

GH (*2H -(-1) 

takes the form R = — " t , and hence it takes large values, diverging at the Rip time t = t s . The function 

f(R), for a large i?, takes the form: 



f(R) ~ R x - a - /2 
2. If 5 - 2\/6 < H Q < 5 + 2\/6, the solution of is given by: 



P{<t>) = («, 



-(h +l)/2 



A COS ( (t 



i In ■ 



(31) 



-Hi + 10H - 1 



B sin (t s — <ft) In ■ 



-tf 2 + 10H - 1 



(32) 



Then, for this choice of the constant Ho, and by means of the equation |(5J| the form of the function f(R) is 
found: 



f(R) ~ i?(^o + D/4 



Aco S r^ln- g ° +1 2 0g0 - 1 



+Bsinfi?- 1 / 2 ln-^- 



10i? - 1 



(33) 



Hence, the expressions (|3"Tj) and (f3"3")l for /(-R) reproduce the behaviour of the Hubble parameter for large t given 
in (|28p . where a phantom accelerated epoch ocurrs, and the Universe ends in a Big Rip singularity for t = t s . As 
is shown, this model is reproduced by (|2T|) for small t when the curvature R is large, and by (|3ip or (|33|) when f 
is large. For a proper choice of the power a, the solution for large t is given by (|31[) . which in combination with 
the solution (|2~T|) for small t, it looks like standard gravity f(R) ~ i? for intermediate On the other hand, for 
negavite powers in (|3Tj) and in combination with ((27)) . this takes a similar form than the model suggested in 0], 
/(-R) ~ R + R 2 + X/R, which is known that passes qualitatively most of the solar system tests. As this is an 
approximated form, it is reasonable to think that this model follows from some non- linear gravity of the sort Rcf.[36], 
which may behave as R 2 for large R. The stability of this kind of models (for a detailed discussion see [13]), whose 
solutions are given by (|27|) and (|33p . is studied in Ref. [34J, where the transition between epochs is well done, and 
then, the viable cosmological evolution may be reproduced by these models. The quantitative study of the transition 
between different cosmological epochs is beyond of the purposes of this work. 

Suming up it has been shown that any cosmology may be reproduced by /(i?)-gravity by using an auxiliary 
scalar field and resolving the equations (jTTJ) and (TT2"|) to reconstruct such function of the Ricci scalar. To fix the free 
parameters in the theory, it would be convenient to contrast the model with the observational data as the supernova 
data by means of the evolution of the scale parameter which is obtained in the models showed above. 



III. F(R)-GRAVITY AND DARK FLUIDS 

In this section the mathematical equivalence between f(R) theories, that could reproduce a given cosmology as it was 
seen above, and the standard cosmology with a dark fluid included whose EoS has inhomogeneous terms that depend 
on the Hubble parameter and its derivatives, is investigated. Let us start with the modified Friedmann equations ((4]) 
written in the following form: 

3H 2 = + V V°/'(i?) - □/'(#)) - 3H , 

-3H 2 - 2H = (uf{R) - l -f{R) - H^j , (34) 
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where we have neglected the contributions of any other kind of matter. If we compare Eqs. (|34j) with the standard 
Friedmann equations (3H 2 = K 2 p and —3H 2 — 2H = n 2 p), we may identify both right sides of Eqs. with the 
energy and pressure densities of a perfect fluid, in such a way that they are given by: 



f'(R) \2 



f(R) + V V°f'(R) - af(R) - 3H 



P = ~2 



f'{R) 



Uf\R)--f{R))-H 



(35) 



Then, Eqs. (|34|) take the form of the usual Friedmann equations, where the parameter of the EoS for this dark fluid 
is defined by: 



w 



P 
P 



f'(R) 



(□/'(#) - l/(iZ)) - H 



■pfe (\f{R) + V V°/'(i?) - Df'(R)) - m 
And the corresponding EoS may be written as follows: 

1 



P 



AH 



f(R) 



V V°/'(i?) 



(36) 



(37) 



The Ricci scalar is a function given by R — 6(2H 2 + H), then f(R) is a function on the Hubble parameter H and 
its derivative H . The inhomogeneus EoS for this dark fluid (|37p takes the form of the kind of dark fluids studying in 
several works (see 50]-[55]), and particularly the form of the EoS for dark fluids investigated in (55| . which is written 
as follows: 

p = —p + g(H,H,H...), where 

g(H,H,H...)=-AH + VoVoOnf'(R)) + (V lrif(R)){Volnf(R)) . (38) 
Then, as constructed in [55j, by combining the Friedmann equations, it yields the following differential equation: 



H 



-g(H,H,H.. 







(39) 



Hence, for a given cosmological model, the function g given in Q38p may be seen as a function of cosmic time t, and 
then by the time-dependence of the Ricci scalar, the function g is rewritten in terms of R. Finnally, the function f(R) 
is recovered by the expression (|38|) . In this sense, Eq. |39|) combining with the expression (l38|) results in: 



dx(t) 
dt 



+ x(t) 2 = H(t) , 



(40) 



where x(t) = WW) . Eq . 

is a type of Riccati equation, that may be solved by a given Hubble parameter. 
Hence, /(-R)-gravity is constructed from standard cosmology where a perfect fluid with an inhomogeneus EoS is 
included. To show this, let us consider the following example: 



H(t) = H t 



H 1 ,H >0. 



(41) 



This model reproduces a Universe that passes through two epochs, a first decelerated phase and a second one accel- 
erated that is identified with the current epoch. To resolve (|40|) for this example, and for simplicity, we study the 
asymptotic behaviour of the model. Then, for small t, the Hubble parameter takes the form H(t) ~ and the form 
of the function f(R) is found by solving Eq. (|4*0|) . to be: 



f(R) 



l 



1 - k/2 



R l-k/2 + x 



(42) 



where k(l — k) = Hi and A is an integration constant. Then, by the function l|42[) the model (|4Tj) for small t is 
reproduced. The analog dark fluid that reproduces this behaviour, may be found by inserting the function (|42p in the 
EoS for the fluid given by l[3"8"j). it yields: 



1 



-4H+1 (fc/2(l + fc/2)|-^ 



(43) 
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where R = 6(2H 2 + H). Then, a perfect fluid with an inhomogeneous EoS given by (|43|) reproduces the asymptotic 
behaviour for small t in an analog description to f(R). Let us now explore the form of f(R) function and the EoS 
of the dark fluid for large t, in this case H(t) ~ Hot, and by resolving the equation (|40|) . the expression for f(R) is 
given by: 



f(R) ~ exp /^oi? 




(44) 



where A is an integration constant. The alternative description in terms of a dark fluid is shown trough its EoS, which 
is calculated as in the case above, by the equation 



1 

-P+ — 



ATT , R " (1 , JE°) - R 

HqR V 4R ; ^/Hji. 



(45) 



which gives an EoS dependent on the Hubble parameter and higher derivatives contained in the Ricci scalar R = 
6(2H 2 + H). Hence, the example considered (|41l) is reproduced in /(i?)- gravity, where its asymptotical behaviour 
is given by the functions calculated above f (|42[) and ([33])). And the same result may be reproduced by a dark fluid 
whose EoS is given by the functions (|43[) and (j4"S"|) . Thus, as it is showed, f(R)-gravity may be written as dark fluids 
with particular dependence on the Hubble parameter and its derivatives through its EoS (|38[) , so the same model may 
be interpretated in several ways. 



IV. DISCUSSIONS 



f(R)-gravity theories may provide an alternative description of the current accelerated epoch of our Universe and 
even on the whole expansion history. As it is pointed in several works (see [39j-(47l| ), one may construct this kind of 
theories in accordance with the local test of gravity and with the observational data, which provide that, at the current 
epoch, the effective parameter of the EoS is close to -1. The next step should be to compare the different cosmological 
tests, as the supernovae luminosity distance or the positions of the CMB peaks, with the F(R) models, but this is 
beyond of the purpose of this paper. In other hand, we have shown two different ways of reconstruct the f(R)-gravity 
in the context of cosmology, in the first one, an auxiliary scalar field is used, and in the second one, the mathematical 
equivalence between f(R)-gravity and dark fluids with inhomogeneus EoS shows that while the expansion history of 
the Universe may be interpreted as a perfect fluid whose EoS has dependence on the cosmological evolution, this effect 
may be caused by the modification of the classical theory of gravity. However, there is not any complementary probe 
to distinguish between both descriptions of the evolution of the Universe, and thus, such kind of modified gravity is 
completly allowed. Hence, f(R)-gravity is an acceptable solution to the cosmological problem, that may provide new 
interesting constraints to look for. 
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